Non-equilibrium current noise in mesoscopic disordered SNS junctions 
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Current noise in superconductor-normal metal-superconductor (SNS) junctions is calculated 
within the scattering theory of multiple Andreev reflections (MAR). It is shown that the noise 
exhibits subharmonic gap singularities at eV — 2A/n, n — 1,2, ... both in single-mode junctions 
with arbitrary transparency D and in multi-mode disordered junctions. The subharmonic structure 
is superimposed with monotonic increase of the effective transferred charge q* = Si(0)/2I with 
decreasing bias voltage. Other features of the noise include a step-like increase of q* in junctions 
with small D, and a divergence Si(0) oc V' 1 ' 2 at small voltages and excess noise Sex — 2el ex , where 
I ex is the excess current, at large voltages, in junctions with diffusive transport. 



The process of multiple Andreev reflections (MAR) in 
Josephson junctions with large electron transparency D 
leads to a very non-trivial behavior of their transport 
properties, most notably, to the subharmonic gap sin- 
gularities in the dc current / at voltages eV — 2A/n, 
n = 1,2, ... , where A is the superconducting gap - see, 
e.g., [1] and references therein. In the past few years, 
experiments with atomic point contacts fabricated by 
the controllable break junction techniqueou have demon- 
strated a remarkable level of agreement between theory 
and experiment in the description of the subgap structure 
in the dc currentem. 

The noise properties of highpjbransparency Josephson 
junctions are also non-triviaB a. The main qualitative 
feature of the noise in the MAR regime is the increase 
of the effective transferred charge S/(0)/2i", where Si(u) 
is the spectral density of current noise, with decreasing 
bias voltage. When the voltage decreases, the number n 
of Andreev reflections in MAR cycles giving the dominant 
contribution to transport increases as 2A/eV '. Since each 
pair of Andreev reflections from the two superconducting 
electrodes transfers the charge 2e across th&iunction, the 
effective noise charge increases as 2A/V.Qm!i] 

Despite the fact that considerable experimental effort 
has been focused recently on measuring noise in diffu- 
sive superpQiuinctor-normal metal-superconductor (SNS) 
structurealj Ej, the noise has been studied theoretically 
only in single- mode ballistic junctions with D = 1, where 
the enhancement of q* manifests itself only at tempera- 
tures T ~ A and disappears at T — > 0. The aim of 
this work was to study the current noise associated with 
MAR in multi-mode junctions with arbitrary mode trans- 
parencies. We find that at finite reflection coefficients the 
large low- voltage noise exists even at zero temperature, it 
exhibits pronounced subgap structure, and is again char- 
acterized by the effective charge on the order of 2A/V. 

In what follows, we use the standard scattering ap- 
proach to the description of MAR in voltage-biased short 
disordered Josephson junctions - see, e.g., [13]. The junc- 
tion is modeled as a normal constriction between two su- 
perconductors with energy gap A. The constriction re- 
gion is assumed to support N propagating electron modes 
and is characterized by the 27V by 2N scattering matrix 
S. Our objective is to find the spectral density Si(ui) of 



current fluctuations 



Sj(«) = / dre WT {I{t)I{t + t) + I(t + r)I(t)) , (1) 



where the bar indicates averaging over the time t, I(t) = 
I(t) — (I(t)) is the deviation of the current from its av- 
erage value, and (...) denotes averaging over the equilib- 
rium distribution of quasiparticles incident on the junc- 
tion. The operator / of the current in the normal con- 
striction is expressed in the standard way through the 
electron operators: 
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Since the superconductors induce electron-hole corre- 
lations in the normal region, the electron operator should 
be written as in the Bogolyubov's transformation: 



*<t = )Xu v b Vj , 



+ siga(a)v*bt ) 



(3) 



The operators 6„. CT describe quasiparticles with spin a in- 
cident on the junction. The quasiparticle index v consists 
of the mode number i, energy e, and the index s = L, R 
indicating whether the quasiparticle is incident from the 
left or right electrode: v = {i, e, s}. Substituting Eq. (||) 
into Q2J) and then into (|l ) we can take two steps in trans- 
forming the expression (fl) for the noise Si(u>). As a first 
step, using the properties of the scattering matrix S (uni- 
tarity and time reversal symmetry) we prove that Si(u>) 
can be presented as a sum of independent contributions 
of different modes. Then, we obtain expression for the 
contribution of a single mode in terms of the amplitudes 
of the forward- and backward-propagating components 
of electron and hole wavefunctions u, v in Eq. (m. The 
four relevant amplitudes are the amplitudes A n , B n of 
Andreev and normal reflection with energy e + 2neV for 
quasiparticle incident at energy e, and the amplitudes 
C„, D n of the normal and Andreev transmission with en- 
ergy e+ (2n + l)eV. All these amplitudes are determined 
by the by-now standard set of recurrence relations that 
describe the process of MAR - see, e.g., [13]. In terms of 
this amplitudes the current noise of a single mode is: 
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Here a^ = a(e + fceF) with a(e) = sign(e)(|e| — (e 2 — 
A 2 ) 1 / 2 )/A being the amplitude of Andreev reflection 
from a superconducting electrode, A n = A n + S n fi/ao, 
and /(e) is the equilibrium Fermi distribution of quasi- 
particles. All quantities with primed indices in Eq. (]4|) 
are functions of e'. 

Despite its apparent complexity, the structure of Eq. 
(Q) is quite simple. The four groups of terms in this equa- 
tion correspond to four different contributions to noise 
coming from the correlations between quasiparticles of 
the same type/different types (e.g., electron-electron or 
electron-hole) incident from the same electrode/opposite 
electrodes. The first group represents the contribution of 
similar quasiparticles coming from the same electrode, 
the second - similar quasiparticles from opposite elec- 
trodes, etc. Results of the numerical calculation of noise 
for zero frequency and temperature from Eq. (4) are 
shown in Fig. 1. 

We see that the main qualitative feature of the noise for 
all transparencies D not too close to and 1 is the subhar- 
monic gap structure at voltages eV = 2A/n, n = 1,2, ... 
that is similar to, but more pronounced than, the sub- 
gap structure in the average current (shown in Fig. 1 
by dashed lines for comparison). Both structures arise 
from the quasiparticle trajectories in energy space which 
hit the two edges of the superconducting energy gap, 
e = ±A, leading to resonant transmission of quasipar- 
ticles through the contact. With the normalization used 
for Sj(0) in Fig. 1 the noise and current curves should 
coincide when the noise is mainly due to uncorrelated 
transport of Cooper pairs. Figure 1 shows that this is 
indeed the case for relatively small transparencies in the 
voltage range A < eV < 2A. At very small transparen- 
cies D — > the noise in the subgap region eV < 2A 
disappears together with the current and the curves ap- 
proach the usual tunncl-Hamiltonian rcsultihr quasipar- 
ticle noise in Josephson tunnel junctionslijHij. This can 
also be checked analytically from Eq. (0) for arbitrary 
frequency u> and temperature T. 
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FIG. 1. Solid lines: Zero-temperature spectral density of 
current fluctuations 57(0) as a function of the bias voltage V 
in single-mode junctions with different electron transparencies 
D. Dashed lines: the average dc current I. 



In the opposite limit of ballistic junction with D = 1 
Eq. (0J reproduces ±he result for ballistic point contacts 
obtained previously!! At oj = 0, T = this result implies 
that, in agreement with Fig. 1, the noise is small for a| 
bias voltages and approaches 16e 2 A/157r?i at eV ^> A 
Figure 1 shows, however, that the curves with finite re- 
flection coefficient R = 1 — D ^ approach the ballistic 
limit R = in a non-trivial manner. As R decreases, 
the noise develops a peak at small voltages. The peak 
becomes narrower and higher with decreasing R and is 
shifted to smaller and smaller voltages, so that at R —* 
it becomes infinitely narrow and disappears. 

This behavior of the noise at small bias voltages 
can be understood quantitatively using the concept of 
Landau-Zener transitions between the two subgap An- 
dreev states developed before for the description of the 
average currentJl3. The two states carry the currents 



±(eA/7i)siny>/2, where ip is the Josephson phase dif- 
ference across the junction, (fi = 2eV/h. The origin 
of the noise in the regime of the Landau-Zener transi- 
tions is the stochastic quantum-mechanical nature of the 
transitions between the two states. In each period of 
the Josephson oscillations the junction either stays on 
one of the Andreev levels and carries the current of the 
same sign during the whole oscillation period, or makes 
the transition between the two states at <p = it so that 
the current changes sign for the second half of the pe- 
riod. The first scenario is realized with the probability 
p = cxp{— irRA/eV} and then the charge 



Qo 



eA 



at sva.ip/2 = — — 
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where To = irh/eV is the period of the Josephson os- 
cillations, is transferred through the junction. The sec- 
ond scenario has probability 1 — p and no net charge is 
transferred in that case. Therefore statistics of charge 
Q transferred through the junction during the large time 
interval T 3> To is characterized by the binomial distri- 
bution with probability p. Definition (yj) then gives: 
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(5) 



The noise (H) reaches maximum at small voltages eV ~ 
irRA and its peak value increases with decreasing R as 
1/R. One can check that Eq. (H) describes accurately 
the low-voltage behavior of the curves with small R in 
Fig. 1. 

Comparing Eq. (J5|) to the expression for the average 
current in this regimetS, I = (2eA/irh)p, we see that 
at small voltages, eV <§; irRA, when p <§; 1 and both 
the average current and noise are suppressed, the effec- 
tive charge q* diverges as l/V, q* = 5/(0) /2J = 2A/V. 
The reason for this divergence is the coherent transfer of 
charge within each cycle of MAR, the number n of An- 
dreev reflections in which is determined by the voltage, 
n ~ 2AjeV . Each cycle transfers the charge en, and 
the current flows through the junction by rare random 
avalanches of Andreev reflections occurring with proba- 
bility p and transferring charge quanta q* = en = 2AJV 
each. Therefore the noise in this regime is the shot noise 
of the large charge quanta q* . 

From the approximate relation q* — en, one can expect 
that q* could exhibit a step-like structure as a function 
of bias voltage due to discreteness of n. In junctions 
with large transparency D ~ 1, however, different cycles 
of MAR contribute to transport simultaneously, and dis- 
creteness of n is washed out. Nevertheless, when D is 
small, the probability of the cycle of n Andreev reflec- 
tions is proportional to D n , so the cycles with smallest 
possible n dominate the transport. In this case the step 
structure in q* is pronounced, as can be seen from Fig. 
2, where the average current and current noise are plot- 
ted as functions of voltage for two values of transparency 



D <C 1 on a logarithmic scale. The noise curves in Fig. 
2 follow closely the current with a shift that increases 
with decreasing voltage and reflects the increase of q*. 
The inset in Fig. 2 shows directly the step structure in 
q* as a function of 2A/V^, and its gradual suppression 
with increasing D. 




FIG. 2. Current noise <Sr(0) and the average current I 
(same as in Fig. 1) on a logarithmic scale for two small values 
of the junction transparency D. Voltage dependence of the 
effective charge q* = Si(0)/2I is shown in the inset. 

In the case of multimode junctions all the results ob- 
tained above should be averaged over the distribution of 
the mode transparencies. Below, we do this for the most 
important case of the SNS junction, where the N region 
is a disordered normal conductor with a large number 
of propagating modes and diffusive electron transport. 
Such a conductor is characterized by quasicontinuous 
Dorokhov's distribution of transparencies (see, e.g., [17] 
and references therein): p(D) — irhG/2e 2 D(l — D) 1 / 2 , 
where G is the normal-state conductance of the N region. 
Averaging the numerical results for noise in single-mode 
junctions with this distribution we obtain the voltage de- 
pendence of the non-equilibrium current noise in the SNS 
junctions, as shown in Fig. 3. We see that the averaging 
preserves the two qualitative features of the noise: sub- 
harmonic gap structure in the subgap region and excess 
noise at eV 3> A. The new feature of 5/(0) is the diver- 
gence at low voltages eV <C A. This divergence arises 
from the contribution of modes with small reflection co- 
efficients and corresponds to the low-voltage noise peaks 
in the single-mode junctions. Averaging Eq. (g|) for these 
peaks with the Dorokhov's distribution we get the ex- 
pression for the low-voltage noise in the SNS junction: 
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(6) 



Equation (||) agrees with Fig. 3 and shows that at 
eV <§; A the noise diverges as 1/y/V. In presence 
of a finite energy relaxation (not accounted for in our 
scattering approach) this divergence wouldi saturate at 
eV ~ tvy, where 7 is the relaxation rate.liS Combining 



Eq. (g) with the similar expression for the dc currently, 
I = G(V A/e) 1 ' 2 , we find that in the diffusive junction 
the effective charge q* has the same 1/V voltage depen- 
dence as in the nearly ballistic single-mode junction (but 
with a different coefficient): 



(I--)-- 
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(7) 



Equation (|7J) represents the main trend in the behavior 
of q* as a function of the bias voltage. The inset in Fig. 
3 shows that on top of this linear growth of q* with 1/V, 
it exhibits weak periodic oscillations associated with the 
subharmonic gap structure. 
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FIG. 3. Current noise Si(0) (solid line) and the average 
current / (dashed line) of a short diffusive SNS junction as 
functions of the bias voltage V at zero temperature. The inset 
shows the effective noise charge q* versus the inverse voltage. 

Next we discuss the excess noise. It can be seen 
from Eq. (4) and easily understood on physical grounds, 
that the leading contribution to the noise at large volt- 
ages eV _S> A coincides with the noise of the normal 
junctionsril, S/(0) = (2e 3 /irh)VD(l - D) in the single- 
mode and S/(0) = 2eGV/3 in the multimode diffu- 
sive junction. In addition to this leading normal term 
that increases with voltage, there is an extra voltage- 
independent contribution S ex to the large-voltage noise. 
This "excess" noise can be seen in Figures 1 and 3 
and is directly related to the process of Andreev re- 
flection. Indeed, one contribution to S ex comes from 
the quasiparticles incident on the superconducting elec- 
trodes at energies e outside the energy gap, | e |> A. 
For such quasiparticles, the probability of Andreev re- 
flection is less than one and the random process of An- 
dreev reflection/transmission introduces extra noise in 
addition to the noise associated with random normal re- 
flection/transmission due to scattering in the normal re- 
gion. Another contribution to S ex is due to quasiparticles 
inside the energy gap, | e |< A. Since each quasiparticle 
that is Andreev reflected in this energy range transfers 



the charge 2e (instead of the "normal" e), Andreev reflec- 
tion amplifies the randomness of the normal scattering 
and increases the noise beyond its normal-state value. 

An expression for the excess noise in a single-mode 
junction as a function of transparency D can be obtained 
from Eq. (|J) using the fact that only MAR processes 
with one Andreev reflection contribute to the noise in 
the eV 3> A limit. Averaging the single- mode expres- 
sion over the Dorokhov's distribution we find the excess 
noise of the diffusive SNS junctions: 



S ex = 2GA(— - 1) = 2el e 



(8) 



The last equality in Eq. (pf) follows from the known ex- 
pression for the excess current of such junctionsEI We see 
that due to a curious numerical coincidence, the excess 
noise and current are related by a regular Schottky-like 
relation, in contrast to the leading normal contribution 
to the noise which is a factor of 1/3 smaller than the 
Schottky expression. 

In summary, we have calculated the non-equilibrium 
current noise in the MAR regime in voltage biased short 
multimode Josephson junctions with arbitrary electron 
transmission properties. The noise exhibits several qual- 
itative features associated with the coherent charge trans- 
fer by MAR cycles: subharmonic gap structure, "quan- 
tization" of the effective charge q*, pronounced zero- 
temperature singularity at low bias voltages, and excess 
noise at large voltages. 
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